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We present a detailed next-to-leading order (NLO) leading twist QCD analysis of deeply virtual 
Compton scattering (DVCS) observables, for several different input scenarios, in the MS scheme. We 
discuss the size of the NLO effects and the behavior of the observables in skewedness C^, momentum 
transfer, t, and photon virtuality, q = —Q . We present results on the amplitude level for un- 
polarized and longitudinally polarized lepton probes, and unpolarized and longitudinally polarized 
proton targets. We make predictions for various asymmetries and for the DVCS cross section and 
compare with the available data. 



I. INTRODUCTION 



In the quest for understanding the structure of hadrons, hard, exclusive lepton-nucleon processes have emerged as 
very promising candidates to further constrain the dynamical degrees of freedom of hadronic matter. Experimentally, 
I such processes are typified by a clear spatial separation of the scattered final state nucleon and the diffractively- 

■ produced, exclusive system, X, i.e. by the presence of a large rapidity gap. The hard scale required for a perturbative 
' analysis is either provided by the spacelike virtuality, = —q^ 3> ^qqb^ exchanged photon, a heavy quark 

■ mass, or by a large momentum transfer to the hadron in the t-channel, t <^ 0. Deeply virtual Compton scattering 
(DVCS) pl-p^, j*{q) + p{P) —^ 7(f?') + p{P'), is the most promising of these processes. One reason for this is 
that on the lepton level it interferes with a competing QED process, known as the Bethe-Heitler (BH) process, in 
which the final state photon is radiated from either the initial or final state lepton. The associated interference 
term offers the unique possibility to directly measure both the imaginary and real parts of QCD amplitudes, via 

' various angular asymmetries. A factorization theorem has been proven for the DVCS process |^,^ which relates the 
experimentally-accessible amplitudes to a new class of fundamental functions, called generalized parton distributions 
O i' (GPDs) p]-p|,p^|-p7[ , which encode detailed information about the partonic structure of hadrons. 
D . GPDs are an extension of the well-known parton distribution functions (PDFs) appearing in inclusive processes such 
as deep inelastic scattering (DIS), or Drell-Yan, and encode additional information about the partonic structure of 
hadrons, above and beyond that of conventional PDFs. They are defined as the Fourier transforms of non-local light- 
cone operators!^ sandwiched between nucleon states of different momenta^, commensurate with a finite momentum 
transfer in the t-channel to the final state proton. As such, the GPDs depend on /oitr variables {X,(^,Q'^,t) rather 
than just two {X, Q^) as is the case for regular PDFs. This allows an extended mapping of the dynamical behavior of a 
nucleon in the two extra variables, skewedness C, and momentum transfer, t. In fact, knowledge of the behavior of the 
GPDs in these two extra variables would allow one to obtain, for the first time, a three dimensional map of the proton 
in terms of its partonic constituents. The GPDs are true two-particle correlation functions, whereas the PDFs are 
effectively only one-particle distributions. They contain, in addition to the usual PDF- type information residing in the 
so-called "DGLAP region" [p^ (for which the momentum fraction variable is larger than the skewedness parameter, 
X > C), supplementary information about the distribution amplitudes of virtual "meson- like" states in the nucleon 
in the so-called "ERBL region" @ (^ < C)■ 



*andreas.freund@physik.uni-regensburg.de 
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^The fact that the produced real photon is an elementary quantum state eliminates the need for further non-perturbative 
information, which is required for example in exclusive vector meson production [Q. This simplifies the theoretical treatment 
considerably. 

■^compared to local operators in inclusive reactions, 
^in inclusive reactions the momenta are the same. 
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A good knowledge of GPDs is required to establish the boundary conditions for a large class of exclusive processes 
calculable in QCD. Unfortunately, reliable perturbative calculations can only be made if t is either small or large, 
confining a sensible comparison between theory and experiment to restricted kinematical regions where either the 
t-dependence is a purely nonperturbative function (small large Q^) as in DVCS or the Q^-dependence is mainly 
non perturbative (small Q^, large t) as in, for example, large-t photoproduction of a real photon (wide angle DVCS) 
pO| . In light of this observation, one might question the practicality of studying and measuring these exclusive 
distributions, given that inclusive PDFs may only be constrained well via a global analysis of a large number of data 
points from numerous experiments. 

It turns out that the GPDs are rather more tightly constrained than one might naively assume . Firstly, they are 
obliged to reproduce the regular PDFs in the forward limit ( ^ ■ Secondly, they are each required to be either 
symmetric or antisymmetric about the point X = in the ERBL region, and they have to obey a polynomiality 
constraint (see for example |^). These are properties which need to be preserved under evolution. Lastly, the 
DVCS amplitudes, and thus certain observables, seem to be very sensitive to the shape of the GPD in the small ( 
region, especially the real part of DVCS amplitudes 24|. Hence, experimental measurements of DVCS observables, 
even of only moderate statistics, appear to give a good opportunity to pin-down the GPDs, given these theoretical 
restrictions. Therefore a careful, thorough and accurate theoretical analysis of DVCS is necessary to understand how 
varying the input GPDs affects the physical observables and the quantitative and qualitative changes in going from 
leading order (LO) to next-to-leading order (NLO) accuracy in perturbation theory. In this paper, we present such 
a NLO analysis, in the MS-scheme, for both polarized and unpolarized scattering, and explore some of the necessary 
issues required to make an optimal extraction of the GPDs from current and future data. 

The physical picture emerging for DVCS is also very interesting in its own right. Several important questions 
immediately arise. What does the energy and Q^-dependence of DVCS reveal about the nature of diffractive exchange 
and how does it compare to other diffractive processes ? Why does DVCS appear to have significant probability in 
the valance region at larger x, i.e. outside of the region in which generic diffraction usually occurs ? Is the physical 
picture for DVCS the same in both regions ? We will attempt a partial answer to these questions in the following. 

This paper gives a comprehensive analysis of DVCS observables at NLO accuracy. Complementary information 
may be found in |^,^,|4|. To bring our analysis right up to date, we introduce unpolarized input GPDs based on two 
of the most recent PDF sets, CTEQ5M || and MRST99 [|6) (in addition to GRV98 (2|], and the older MRSA' |2|] 
used in our earlier publications). This allows us to push our input scale for skewed evolution down to Qo = 1 GeV. 
For the purposes of comparison we also present the results for DVCS observables obtained using our earlier input 
models. Various computer codes used in our analysis our available from the HEPDATA website 

This paper is structured as follows. In section || we reiterate the kinematics of DVCS and BH and define the 
DVCS observables, i.e. various measurable angular asymmetries and cross sections. In section III we describe the 
various input GPDs and present numerical results for the NLO evolution of the new sets. Section fv\ discusses 

how the various DVCS amplitudes are produced via convolution integrals of GPDs with coefficient functions. We 
present the and ^-dependence of the unpolarized amplitudes for the new sets graphically. In section |^ we give 
our predictions for the DVCS observables in C, and t, for various polarizations of probe and target, as well as 
discussing their implications. We compare our results with the cur rently available data and with other theoretical 
predictions in section VI. Finally, we briefly conclude in section 



VII. 






a) b) c) 

FIG. 1. a) DVCS graph, b) BH with photon from final state lepton and c) with photon from initial state lepton. 



II. KINEMATICS AND OBSERVABLES FOR DVCS AND BH 
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A. Kinematics and frame definition 



The lepton level process, e'^{k, n) N{P, S) e*(A:', k') N{P', S') j{q', e'), receives contributions from each of the 
graphs shown in Fig. |^. The corresponding differential cross section is given byQ: 

where the square of the amplitude receives contributions from pure DVCS (Fig. la), from pure BH (Figs, lb, Ic) and 
from their interference (with a sign governed by lepton charge), 

l^^l^ = I]K',S',e' \TdVCS? + {'^DVCS'^BH +TdVCs'^Bh) + \'^Bh\'^ ■ (2) 

The DVCS amplitude is given by 



-^Dvcs = ±p;T^''uik'hMk) { ^ t , (3) 



e" 



where q ~ k — k' , e* is the polarisation vector of the outgoing real photon, and the hadronic tensor, T^i/, is defined 
by a time ordered product of two electromagnetic currents: 

T^.(g- F, P')=^J dxe'^-HP', S'\Tj,ix/2)j,i-x/2)\P, S) , (4) 

where q = {q + q')/^. This hadronic tensor contains twelve]^ independent kinematical structures for a spin-1/2 
target. Here, we restrict ourselves to the twist-2 part of r^,y and drop all other terms of either kinematical or dynamical 
higher twist. From the structure of the operator product expansion (OPE) one can immediately conclude that, to 
twist-2 accuracy^, one has the following form factor decomposition: 

T^y (9, P, A) = -5^^ ^p-^ - ie^,yp„ -^p—z , (5) 

where P = (P + P')/2 and the gauge invariant tensors g^^i, — V^pgpaVau and e^uaf} = V^pepaapVau are constructed 
through the projection tensor Vpv = g^v — qtiq'v/l ' 1' ■ The vector V^^ and axial-vector are expressed, again to 
twist-2 accuracy, through the following form factor decomposition, 

V,, = u{p', s') {n,^, - £1'-^^) uiP, s) , (6) 

A,p = U{P', S') (^Hi7p75 - ^1 ^) U{P, S) , (7) 

where U, U are spinors for the incoming and outgoing hadron state, A = P — P' is the momentum transfered from 
the hadron^ and M is the hadron mass. The various Lorentz structures have associated amplitudes: TCi,£i are 
unpolarized helicity non-flip and helicity flip amplitudes, respectively, and Ti.i,£i are their polarized counterparts. 
These amplitudes are expressed, via the DVCS factorization theorem as convolutions of a hard scattering 

coefficient function and a GPD. They depend on the following Lorentz-invariant variables: 

C = J^, Q' = ~f, ^ = A2^(P-P')^ 

2P ■ q 



■^In this section we follow closely the notation of 

®12 = ^ X 3 (virtual photon) x 2 (final photon) x 2 (initial nucleon) x 2 (final nucleon). The reduction factor 1/2 is a result 
of parity invariance. 

®We drop a twist-2 contribution arising from a double helicity flip of the photon, i.e. going from helicity -1-1 to helicity — 1 
or vice versa, which is suppressed in Us since this double flip can only be mediated by gluons. Thus when we speak of twist-2 
contributions we really mean twist-2 modulo this double flip contribution. 

'^Note that there is a relative minus sign between our definition of A and that of [KlJ (our definition is the same as r in p2]). 
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which are related to the experimentally accessible variables, C = ^bj = —q'^/{2P ■ q) and = — g^, used throughout 
this paper, via 



The BH amplitude is purely real and is given by the sum of the graphs in Fig. lb and Fig. Ic: 

TBH = -je';L^''J,, (9) 

with the leptonic tensor 

L^, - uik', k') [7^(^ + A)"'7. + 7.(r - A)"'7/.] n{k, n) , (10) 

and the hadronic current 

.h - U{P',S') (^Fi(i)7, - iF2{t)a,r^^ U{P,S) , (11) 

where Fi, F2 are the Dirac and Pauli form factors, respectively, normalized such that Ff{0) — 1, -Ff (0) = — —1.79, 
and F"(0) = 0, i^^(O) = k„ = —1.91 for proton, p, and neutron, n. These are known from low energy exclusive 
scattering and have been parametrised, using dipole formulas for small t, as linear combinations of the electric and 
magnetic form factors: 

^1= . , t ' ^2[t) = . , t ' (1^) 

with 

= ^ = ^ = , ' Gl;(t) = 0, (13) 
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and my = 0.84 GeV ||. 

Following [Q, we chose to work in the 7*P frame with the proton at rest such that the direction of the four 
vector q = k — k' (i.e. the virtual photon direction for the DVCS graph) defines the negative z-axis ^. Without 
loss of generality we can choose the incoming electron to only have a non-zero component along the positive x- 
axis in the transverse (x — y) plane: k = (feg, fco sin^e, 0, /cq cos^e), q = {qQ,0, —\q3\), P = (M,0,0,0) and P' = 
(Pqi l-F"l cos0sinf?/f , |P'| sin(/)sin0//, COS0//), where (j) is the azimuthal angle between the lepton (a; — z) and 
hadron scattering planes. A dependence on the angle (j) can be either induced kinematically, or through the phase of 
hadronic amplitudes as will be shown below. The spin vectors of the nucleon target for longitudinal and transverse 
polarizations are given by 

± _ A / Q' 



S^^ ===[qo--^,0,0,-\q3\)] , = (O,cos$,sin$,0), (14) 



where A = ±1. Note that t = (P — P')^ has a kinematical minimum value which is given by 

M'^xh ]\/P 
Wr = -^--^(l + (^(^)). (15) 

The motivation for using this frame is that the frame-dependent expressions for the s- and it-channel BH propagators 
appearing in eq.(p^ (cf. Figs, lb, Ic, respectively) have a particularly simple Fourier expansion^ in the angle 



This frame is related to the center-of-mass system of ^ by a boost of the hadron in the z-direction. 
^In a different frame is used which induces a more complicated ^-dependence than one finds in our frame. 
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(j). The dimensionless forms Pi = = *^q2^ and V2 = q2 = 1 ^ are (up to corrections of order 

0{{-t/Q%{M^/Q^))): 

_{l + 2^ '°"° -" cos^- ^(1 - , - x(2 - ,)) - 2— (2 - ,)} , 

. - ,) + 2^ "-'°-'"^;-"'-'' Ls^- - x,2 - „)) - 2ilff!,2 - .)} ^ (M) 

The product V1P2 appears in the BH and interference expressions and thus induces an additional kinematical (rather 
than hadronic) ^-dependence. In certain kinematical regions, this additional ^-dependence can fake certain hadronic 
(^-dependences |p4|| , it can also lead to unwanted contributions in certain 0-asymmetries, as discussed below. One way 
out of this dilemma is to weight DVCS observables with Pi 7^2, leaving only the pure hadronic 0-dependence (exploiting 
the orthogonality of cos (/) and costo'0 for integer m! ^ 1). We will not do this here because such a weighting requires 
a good <j) resolution not available for the present data and also because the pure twist-2 contributions may well explain 
most of the observed data without needing to take twist-3 or higher contributions into account. Nevertheless, studying 
weighted DVCS observables should be done as soon as good experimental </)-resolution is available. 



B. DVCS observables: differential cross section and asymmetries 



After performing the phase space integration in eq.(0), the triple differential cross section on the lepton level is 
given by 



dxbjdQ'^d\t\ Jq dxbjdQ^d\t\d(j> SttQ'* 

The twist two expressions for the DVCS squared, interference and BH squared terms, for all probe and target 
polarizations, required for eq. (p7|) are very similar to eqs. (24-32) of but with the full expressions for the BH 
propagators of eq.(16) included to reinstate the correct y- and 0-dependence (a correction factor of —{l — y)/{y^'PiP2) 
should be applied to eqs. (27-32) of ||l|). 

From the pure DVCS piece (following the usual single photon exchange flux factor convention adopted in p5[), 
changing variable from xtj to y and integrating over t, one may define the virtual-photon proton cross section via 

dcr^'He^P ^ e^jp) ^ a.^. (1 + (1 - y)^) . ^ p. .^g. 
dydQ^ 27TQ^y ^ ' ' ^ ^' 

We give predictions for a{'y*P 7P) and compare with the recent experimental data from the HI CoUab. in 
section 
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We will now define various DVCS observables, in terms of a list of asymmetries in the azimuthal angle 0: 

• The (unpolarized) azimuthal angle asymmetry (AAA) , measured in the scattering of an unpolarized probe on 
an unpolarized target, is defined by 

AAA = , (19) 

/o " d^daDVCS+BH 



where da^^ is the pure BH term. 

The single spin asymmetry (SSA), measured in the scattering of a longitudinally polarized probe on an unpo- 
larized target, is defined by 

SSA = , (20) 

^^"^ d(j){dcjDVCS+BHA j^^„DVCS+BH,l) 

where Act = da^ —da^ and t and [ signify that the lepton is polarized along or against its direction, respectively. 
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The asymmetry of an unpolarized probe on a longitudinally polarized target (UPLT) is given by: 

where Actlt = rfo"! — da^ with f and | signifying that the target is polarized along or against the +z-direction, 
respectively, corresponding to A = =pl in eq.(n4]). 



The asymmetry of an unpolarized probe on a transversely polarized target (UPTT) is given by: 



UPTT = i-^ ^ , (22) 

J^"" d<i){daEyc:s+BH + ^^DVCS+BH-^ 

where Acttt — da—, — da^ with — > and ^ signify that the target transverse polarization vector, S'^ , points 
along the +x and —x directions (i.e. $ = 0, tt), respectively. 

The charge asymmetry (CA) in the scattering of an unpolarized probe on an unpolarized target: 



CA = , (23) 

/o'" dcf>id+aD^CS+BH + d-^DVCS+BH-^ 

where Ad^' a — d^a — d^ a corresponds to the difference of the scattering with a positron probe and an electron 
probe. 

• The charge asymmetry with a double spin flip of a longitudinally polarized probe on a longitudinally polarized 
target (CADSFL): 

#Ad^af^^^+^^'^^ - d0Ad^af^^^+^^'^^ 

CADSFL = — , (24) 

" d</.(d+af ^^^+^^ + d-af ^^^+^^) 

where Ad'-^af;!? — d+cr| — d^crj — IS.d^ a with d+cr| corresponding to a positron beam polarized along its own 

direction scattering with a target with its polarization vector having a positive z-component, d~ corresponds 
to an electron beam polarized against its own direction scattering on target with its polarization vector having 
a negative z-component and with Ad*^cr having the same meaning as for CA. 

• The charge asymmetry with a double spin flip of a longitudinally polarized probe on a transversally polarized 
target (CADSFT): 

C')^ d0AdCa?^^^+^^'^^ - f-'^ d0Ad^a?^^^+^^'^^ 
CADSFT = — , (25) 

Jo'" d4,{d+a^VCS+BH + d-„DVCS+BH^ 

where Ac?'^ct|^|? — rf+al^ — d^ai_ — Ad^a, with d+cl^ corresponding to a positron beam polarized along its own 
direction scattering on a target with a polarization vector pointing in the +x direction, d^a^ corresponds to 
an electron beam polarized against its own direction scattering on a target with a polarization vector pointing 
in the —x direction, and Ad'^cr having the same meaning as in CA. 

The definitions above make the asymmetries directly proportional to the real part of a combination of DVCS 
amplitudes, in the case of the AAA, CA, CADSFT and CADSFL, and to the imaginary part of a combination of 
DVCS amphtudes for SSA, UPTT and UPLT. If one forms the proper combinations from eqs. (24-32) of jSltl, with 
the correction factor included, one observes that for small- a; DVCS observables the information from transversally 
polarized targets is redundant to the information from longitudinally polarized targets as far as the information on 
the real and imaginary part of DVCS amplitudes is concerned. For large x, this is, strictly speaking, no longer 
true! However, higher twist corrections, especially in the normalization of the asymmetries, will make extraction of 
information on individual amplitudes virtually impossible. For this reason we will focus only on DVCS observables 
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which may be obtained using an unpolarized or longitudinally polarized target. Note that for small x and t, these 
combinations of amplitudes reduce to just the unpolarized (for AAA, CA and SSA) or polarized (for UPLT and 
C ADSFL) helicity non-flip amplitudes . Note also that the definition of the AAA is different from the usual one 
(see e.g. the first reference of [^) and is designed to ensure that the numerator contains only the interference term 
and is thus directly proportional to a DVCS amplitude. This slight change was necessary since, on inspection, it was 
realised that the pure BH contribution to the numerator does not vanish when the <j> integrations in the numerator of 
eq.(|l^) are carried out (due to the correction factor, oc I/V1P2, applied to eq.(27) of Hence the BH contribution 

needs to be subtracted from the differential cross section in order to have an asymmetry which is directly proportional 
to the real part of a hadronic DVCS amplitude. 

Note that in the following we will always assume a positron probe, except in the case of charge asymmetries where 
one needs both positron and electron. Thus for the corresponding electron observables the overall sign in the results 
we will quote below has to be reversed. 

At this point, we wish to make a few general comments on the (^-behavior of the lepton level expressions given in 
eqs. (30-32) of (modified by the propagator factors). Firstly, one can cleanly separate the real and imaginary parts 
of DVCS amplitudes using their different (j) behavior either by taking moments with respect to a function in (j) (usually 
either sine or cosine) thereby projecting out the unwanted contributions, or, equivalently, by forming asymmetries in 
the angle 0, as we did above. Secondly, as was observed in the (/)-dependence of the expressions also allows a 

clean separation into different twist contributions. For example (see [0), the real and imaginary parts of twist-2 can 
be separated from one another by using (/)-moments, since the real and imaginary parts of twist-3 amplitudes have 
a (/)-dependence which is very different from that of twist-2: e.g. cosc/) (twist-2) and cos2(/) (twist-3) Thus 
DVCS allows the real and imaginary parts of hadronic twist-2 and twist-3 amplitudes to be isolated for the first time 
within the same experiment, simply by using different moments or asymmetries in the azimuthal angle 0. Therefore, 
experimental upgrades which enhance the instrumentation in the forward region (see e.g. for a HI Collaboration 
proposal) are vital to maximise the physics scope of these experiments. We note the asymmetries as defined above 
do not necessarily require a good resolution in the angle (/», it is sufficient to specify which hemisphere in (/> a given 
event corresponds to. A reasonable sample of events will then be sufficient to measure the asymmetries. If one wishes 
to produce asymmetries weighted for example with a sine or cosine, then a good phi resolution is of course required. 
This concludes our comments on DVCS/BH kinematics and DVCS observables. 



III. GPDS AND INPUT MODELS 
A. Symmetries and representations of GPDs 

For the definition of our input GPDs we follow precisely the prescription given in pl] |. GPDs result from matrix 
elements for quark and gluon correlators of unequal momentum nucleon states and may be defined in a number of 
ways. Following |^,^, we initially chose a definition which treated the initial and final state nucleon momentum 
(P, P' , respectively) symmetrically by involving parton light-cone fractions with respect to the momentum transfer, 
A = P — P', and the average momentum, P = (P + P')/2. The inherent symmetries of the matrix elements are clearly 
manifest in associated symmetries of the GPDs. We then shifted to a definition js?) based on light-cone fractions 
of the incoming hadron, for the purposes of evolution and a direct comparison with conventional PDFs and with 
experiment. We discussed the manifestation of the symmetries in this representation and explicitly illustrated their 
preservation under evolution. 

Matrix elements of the non-local operators are defined on the light cone and involve a light-like vector (z^ = 0). 
They can be most generally represented by a double spectral representation with respect to P ■ z and A • z [|l|,||j22j 
(see eq.(4) of [0). In accordance with the associated Lorentz structures, the non-singlet, singlet and gluon matrix 
elements involve functions corresponding to proton helicity conservation (labelled with F) and to proton helicity flip 
(labelled with K) which are collectively known as double distributions. Henceforth, for brevity, we shall only discuss 
the helicity non-flip parts explicitly. However, the helicity flip case is exactly analogous. The D-terms in eq.(4) of pl| ] 
correspond to resonance-like exchange p^js^ and permit non-zero values for the singlet and gluon matrix elements 
in the limit P • z — > and A • z / 0, which is allowed by their evenness in _P • z. 

By making a particular choice of the light-cone vector, z'', as a light-ray vector (so that in light-cone variables, 
z± = zo±Z3, only its minus component is non-zero: z'^ = (0, z_, 0)) one may reduce the double spectral representation 
of cq.(4) of defined on the entire light-cone, to a one dimensional spectral representation, defined along a light 
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ray, depending on the skewedness parameter, ^, defined by 



C = A-z/2P-z = A+/2F+, (26) 

which is equivalent to our definition in Sec. || (cf. eq.(||)). The resultant GPDs are the off-forward parton distribution 
fmictions (OFPDFs) introduced in 



H{v,i,t)^ f dx' [ ^ ^ dy'5{x' +iy' -v)F{x\v\t), (27) 

J ~l J -l+\x'\ 

where v e [—1, 1]. In terms of individual flavor decompositions, the singlet, non-singiet and gluon distributions are 
given through 

a 

H''{v,O^H%v,0±HS{-v,0, (28) 

where the upper (lower) signs corresponds to the unpolarized (polarized) case. Note that the symmetries which hold 
for the matrix elements can change for the H^s, due to the influence of the P ■ z, z factors in eq.(4) of In 
particular, the unpolarized quark singlet is antisymmetric about w = 0, as are both D-terms, whereas the unpolarized 
quark non-singlet and the gluon are symmetric. The opposite symmetries hold for the polarized distributions. The 
helicity flip GPDs, are found analogously (double integrals with respect to the Ks) and similar reasoning establishes 
their symmetry properties. 

As in pT[ | , we shall make the usual assumption that the t-dependence of all of these functions factorizes into implicit 
form factors. One should bear in mind that in order to make predictions for physical amplitudes (for t ^ Q) these 
form factors must be specified. Note that the assumption of a factorized ^-dependence, as a general statement, must 
be justified within the kinematic regime concerned. It appears to be valid at small x and small t, from the HERA 
data on a variety of diffractive measurements. However, it appears not to hold for moderate to large t and larger x 

For the purposes of comparing to experiment it is natural to define GPDs in terms of momentum fractions, X G [0, 1] , 
of the incoming proton momentum, P, carried by the outgoing parton. To this end we adapt the notation and 
definitions of introducing two non-diagonal parton distribution functions (NDPDFs), JF* and J^^, for flavor, a: 

where Vi G [— ^, 1], V2 ^ [—1, £] (see Fig. 4 of [0), C = A+/P+ is the skewedness defined on the domain C, G [0, 1] such 
that C ~ C/(2 — C) and C = x^j for DVCS (this definition is equivalent and the relations are the same as those in Sec. 
H). The transformations between the Wi,f2 and Xi^X2 are given implicitly in eq. (p9|) , the inverse transformations 



are: 



Xi - C/2 C/2 - X2 

For the gluon one may use either transformation, e.g. 

rnx.o-"^^. (31, 

There are two distinct kinematic regions for the GPDs, with different physical interpretations. In the DGLAP |T^ ] 
region, X > (\v\ > ^), !F'^{X,C,) and J-^{XX) are independent functions, corresponding to quark or anti-quark 
fields leaving the nuclcon with momentum fraction X and returning with positive momentum fraction X — C^. As such 
they correspond to a generalization of regular DGLAP PDFs (which have equal outgoing and returning fractions). 
In the ERBL region. A" < C (|u| < ^), both quark and anti-quark carry positive momentum fractions (A, C, ~ X) 
away from the nucleon in a meson-like configuration, and the GPDs behave like ERBL [ p^ distributional amplitudes 
characterising mesons. This implies that J-"^ and J-^ are not independent in the ERBL region and indeed a symmetry 
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is observed: ^'(C — X,() = (which directly reflects the symmetry of H''{v,£,) about v = 0). Similarly, the 

gluon distribution, JTf, is DGLAP-like for X > C and ERBL-like for X < (. This leads to unpolarized non-singlet, 
■pNS,a _ j:q,a _ j7q.a ^ ^^j-^j gluon GPDs which are symmetric, and a singlet quark distribution J^^ ~ X^a-^''''' + -^^'^ 
which is antisymmetric, about the point X ~ C,/2 in the ERBL region. Again the opposite symmetries hold for the 
polarized distributions. 



B. GPD input models 

For our input models we we follow precisely the procedure given in section III of which is based on Radyushkin's 
ansatz p(][| for GPDs. Here we briefly describe some salient features which are required for the discussion and give 



various technical details not included in |21|. The input distributions, J-'^'^'^ {XX, Qo), at the input scale, Qq, have 
the correct symmetries and properties and are built from conventional PDFs in the DGLAP region, for both the 
unpolarized and polarized cases. These input NDPDFs then serve as the boundary conditions for our numerical 
evolution. 

Factoring out the overall t-dependence we have the following integral relations between the double distributions 
and NDPDFs for the quark and antiquark: 



l-C/2 i-i J-,+|..| ' ' ^" " (l-C/2) ' 
C) ^ ^ /' 4.' 4y'^ + iy' - ^^0^ - (32) 

with 1 > ui > — ^, —1 < V2 < S, and a similar relation for the gluon (for which one can use either vi or U2). 

Following |3|,^ we employ a factorized ansatz for the double distribution where they are given by a product of a 
proflle function, tt*, and a conventional PDF, {i = q,g). The profile functions are chosen to guarantee the correct 
symmetry properties in the ERBL region and their normalization is specified by demanding that the conventional 
distributions are reproduced in the forward Hmit: e.g. T^{X,( ^ 0) ^ f^{X). The exact i-dependence will be 



specified in Sec. IV, since it depends on whether one is dealing with hclicity-flip or helicity non-flip amplitudes, 
unpolarized or polarized in origin. 

In we specified two particular forward input distributions for the GPDs by using two consistent sets of inclusive 
unpolarized and polarized PDFs, i.e. GRV98 ^ and GRSVOOj^ ^ with A^l^o = 246 MeV, and MRSA' |2|] 
and GS(gluon 'A') || with A^'^^d^^' = 231 MeV at the common input scale Qq = 4 GeV^ and A'^qco'^^ = 174 MeV 
for both sets. These pairs of unpolarized and polarized sets were consistent in the sense that the unpolarized PDFs 
were used to constrain the respective polarized PDFs and both use the same choices for Aq^d, etc in the evolution. 

In order to bring our analysis more up-to-date|^ and to further investigate the input model dependence, we relax 
our (rather weak) consistency requirement in this paper and use two additional contemporary M5'-scheme unpolarized 
sets. We use CTEQ5M [H and MRST99 ||§ in conjunction with the GRSVOOR polarized set at the common input 
scale of Q§ = 1 GeV^ (with A^fcD^^^ = 326 MeV for CTEQ5M0, and A^^i^^^ = 300 MeV for MRST99). For the 
LO evolution of these sets, we have for both models A'^^^^'^^ = 192 MeV. 

Having defined this particular input model for the double distribution one may then perform the y'-integration in 
eq.(|3|) using the delta function. This modifies the limits on the x' integration according to the region concerned: for 
the DGLAP region X > ( one has: 

^9^-(X,C) = 7 / dx'TT-J { x',^-^—^] q"-{x') . (33) 



C 7^ 



the "standard" scenario with an unbroken flavor sea. 
^'^In particular the MRSA'/GS(A) set from 1995 is based on rather old data, 
^■^the "valence" or broken flavor sea scenario. 

'^^we use the FORTRAN code supplied by Pumplin at the input scale of Qo ~ 1 GeV, to feed into our double distribution 
code. Despite the fact that this code is not recommended for use at such a low scale we found that the results matched very 
smoothly onto the results at higher scales. 
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For the anti-quark, since V2 = — wi one may use eqs.( p9| , |3^ ) with V2 —vi, and, exploiting the fact that f^x) 
—q{\x\) for cc < 0, one arrives at 



C 



X , ■ 



-Vl — X 



r(lx'l). 



(34) 



In the ERBL region {X < (, \v\ < £,) integration over y' leads to: 



r „zi±i 



dx'ir'^ x' 



Vl — X 



dx'-K'i X' 



1+5 



Vl — X 



q\\x'\ 



(35) 







(36) 



The non-singlet (valence) and singlet quark combinations are given by: 



^ _ . (38) 

We implemented eqs.(^,^) using eqs. ( |33| , p^ ) and eqs.( |35|j3^ ) for the DGLAP and ERBL regions respectively, 
employing an adaptive Gaussian numerical integration routine on a non-equidistant grid |^. We comment further on 
our usage of grids and integration routines below. 

The integration ranges in eqs.( ^ , ^4[ ) and eqs.( |35| , |36| ) sample the input PDFs all the way down to zero in x' . 
Generally speaking, the providers of PDF sets issue programs that only allow their PDFs to be called for x' greater 
than some minimum value. This is partly for technical reasons but also partly because the PDFs have not yet been 
well constrained by inclusive data in the very small x region, which corresponds to very high centre-of-mass energies. 
For the implementation of MRST, MRSA' and GS we are fortunate to have access to analytic forms themselves at 
the input scale^ and we simply extrapolate these into the very small x region. For GRV98 and GRSVOO this is not 
the case and it was necessary to perform fits, at the Qg-scale concerned, to the small x behavior of these sets for 
values of x where they are available and then extrapolate these fits into the very small x regime^. In investigating 
this issue we noticed that if PDFs (and the quarks in particular) were too singular the integrals for the individual q 
and q were divergent (although this divergence is cancelled in forming the singlet, some regulation of the limit a;' — > 
was required). This issue is not merely of technical interest. The physics message is clear: the GPDs defined above, 
using Radyushkin's ansatz and hence the DVCS observables are very sensitive to the behavior of the PDFs in 
the small x region and also to their extrapolation to extremely small x. To the best of our knowledge this has not 
been explicitly pointed out before. 

The unpolarized singlet also includes the D-term on the right hand side of eq.(^) (in principle there is also an 
analogous term in the unpolarized gluon {D'^ in eq.(4) of |^l[]), but we choose to set it to zero, since nothing is 
known for the gluon D-term except its symmetry.). We adapt the model introduced in p8| for the unpolarized singlet 
D-term, which is based on the chiral-quark-soliton model. This D-term is antisymmetric in its argument, i. e. about 
the point X — C,/2 (in keeping with the anti-symmetry of J-^ ^ and about u = 0). It is non-zero only in the ERBL 
region and hence vanishes entirely in the forward limit. In practice it only assumes numerical significance for large C 
(see Fig. 6 of H). 



This uses up most of the computing time for an evolution run. 
^^ior CTEQ5M we use Pumplin's code || 

^^Having tried several forms we eventually settled on fits of the type f{x) — f{xi) {x/xiY (1 & log(x/a::i) + c (x — a;i)), 
with xi being the minimum value of x available and with the power a constrained to be greater than zero to allow convergence 
as a; ^ 0. 
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C. GPD evolution 



The input GPDs must now be evolved in Q^, using renormahzation group equations, in order to make predictions 
for DVCS ampUtudes at evolved scales. The input GPDs, as previously defined, are continuous functions which span 
the DGLAP and ERBL regions, and evolve in scale appropriately according to generalised versions of the DGLAP or 
ERBL evolution equations. Note that the evolution in the ERBL region depends on the DGLAP region, i.e. there is 
a convolution integral in the ERBL equations spanning the DGLAP region [(, 1], whereas the DGLAP evolution is 
independent of the ERBL region. As the scale increases, partons are pushed from the DGLAP into the ERBL region 
simply through momentum degradation, but not vice versa. The ERBL region thus acts as a sink for the partons. 
Hence, in the asymptotic limit of infinite Q^, we recover the simple asymptotic pion-type distribution amplitude in the 
ERBL region and a completely empty DGLAP region. This will have strong implications for the DVCS amplitudes 
and cross section. Note that in performing the evolution we assumed that the i-dependences of the quarks and giuons 
(which mix under evolution) are the same and factorize such that they do not influence the degree of mixing under 
evolution. Otherwise, any assumed t-dependence will be modified by the QCD evolution, complicating calculations. 
In fact the ^-dependence of quarks and giuons should be different, however, since we study DVCS only at small t, the 
differences in their ^-dependence should be small. This unresolved problem of t-dependence mixing will be addressed 
in another paper. 

The renormahzation group equations (RGEs), or evolution equations, for the DGLAP region and ERBL regions, 
involve convolutions of GPDs with generalised kernels. They are implemented in a FORTRAN numerical evolution 
code. In the DGLAP region the quark flavor singlet and the gluon distributions mix under evolution according to 
generalized DGLAP kernels^^ taken from The flavor non-singlet (NS) quark combinations do not mix under 

evolution. Note that in order to do the full evolution and afterwards extract the various quark species separately one 
needs to solve two separate evolution equations. One for a symmetric combination q"^ — q"" + q'^ — 1/Nf X]a(9° + 
and one for an antisymmetric combination g° = g° — A single quark species, i.e. quark or anti-quark in the 
DGLAP region, or just a singlet or non-singlet quark combination in the ERBL region, can be extracted the following 
way: 

=i(g^±g^+/), with qS = 1-J2(^q- + q-) ^ (39) 



in the DGLAP region and 



qS'- = (ql+qS)^ and g^^'"^ = , (40) 



in the ERBL region. This procedure was adopted in our FORTRAN code. 

A numerical implementation of the convolution integrals of the RNG equations involves specifying a treatment of 
the integrable endpoint singularities. This is achieved via the following definition of the +-distributions (in this case 
we chose to apply it to the "whole kernel" ) in the DGLAP region: 



(41) 



and accordingly implemented in our code. Note that the lower limit of the first integral in the last bracket is C/Uj 
which is not necessarily a grid point. Since we initially used an equidistant grid in the integration variable]^ we 



^'^Note that in the second reference of there was a typographical error in Eq. (188) where the overall sign of the polarized 
pure singlet term in the QQ sector should be — so as to be consistent with Eq. (178). In another typographical error, the 
factor 3 in the first term of the square bracket in the second line of Eq. (194) of the same reference, i.e. the equation for the 
unpolarized GQ kernel, should be replaced by 3(. These mistakes were properly corrected in the implementation of the kernels 
in the GPD evolution code. 

^*Our method of integration is the following: we first introduce an equidistant grid which we then stretch both in the ERBL 
and DGLAP regions with particular transformation functions in order to be able to treat the important region around ^ and 
more accurately. We then compute the Jacobian of this transformation for the inverse transformation we need. On the 
non-equidistant grid we compute the input distributions and then the kernels. Using the Jacobian we transform back onto the 
equidistant grid and perform the convolution integrals using an equidistant grid integration routine like a semi-open Simpson 
to account for the remaining integrable singularities at y and i^. 
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needed to resort to a slower integration routine for a non-equidistant grid such as an adaptive Gaussian integration 
routine. 

In the ERBL region, the quark singlet and gluon again mix under evolution. The generalized ERBL kernels may 
also be found in The +-distribution, again applied to the whole kernel, takes the following form in the ERBL 
region: 



dzF(|,y^^(z,C) = 



dzv (^|,^ ) [T{z,c)-HC,0]+HC,C) 

1 

dzV 



+ 



C 



y z 



rf^V^ 



c 



y z 
C'C 



dzV 



y z 

C'C 



y 



dzV 



y z 

C'C 



[.F(z,C)-^(C,C)]+^(C,C) dz[V 



c c 



y z 

C'C 



z y 



V 



C C 



z y 

C'C 



(42) 



where the terms have been arranged in such away that all divergences explicitly cancel in each term separately and 
only integrable divergences, as in the case of forward evolution at the point y = xtj, remain. The bar notation means, 
for example, z/ C, — 1 — z / C,. In the code, the integrations were dealt with in an analogous way to the DGLAP region. 
Note that for the solution of the differential equation in Q^, we adopted the CTEQ-routines, which are based on a 
Runge-Kutta predictor-corrector algorithm. 




X 



0.0 0.1 0.2 0.3 0.4 0.5 
X 



FIG. 2. Unpolarized NLO input and evolved singlet quark and gluon GPDs at small and large skewedness. The solid curves 
are the input MRST99 GPDs at Qo = 1 GeV, the dotted ones show them evolved to Q = 10 GeV. The dashed curves are input 
CTEQ5M GPDs at Qo = 1 GeV and the dashed-dotted ones show them evolved to Q = 10 GeV. The quark singlet is scaled 
by a factor of 10~* at = 0.0001 and by 10~^ &t — 0.1. For C, — 0.1 the symmetry of the gluon GPD and the anti-symmetry 
of the singlet quark GPD are apparent about the point X — C,/2 = 0.05. 

In Fig.(||) we plot the input distributions for the "new" unpolarized inputs (CTEQ5M and MRST99) at their input 
scale, and at an evolved scale to demonstrate the evolution effects in NLO. 



IV. DVCS AMPLITUDES 



A. Convolution formulae 



In |24|] we investigated unpolarised and polarised DVCS amplitudes in detail. For convenience in the following 
equations we introduce the notation V for the unpolarized case and A for the polarized case, where applicable, they 
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take the upper and lower signs, respectively]^ The factorization theorem for DVCS proves that the amplitude 
takes the following factorized form in the non-diagonal representation (up to power-suppressed corrections of 0{1/Q)): 



2-C 



v c 



1 /2-C 



Nf V C 



c 

2 _ „1 



1 + ^^('^)^^M(X,C,M^t) 



2X Q2 

TV 



(43) 



Note that the second integral is purely real and does not need a -fie prescription since there is no divergence of the 
coefficient function in the integration interval. Also note the opposite sign structure in the quark singlet and gluon case 
due to opposite symmetries of the quark singlet and gluon. For our numerical calculations, we set the factorization 
scale, fi^, equal to the photon virtuality, (in we studied the effects of its variation and found them to be 
rather mild; the associated uncertainties are less than those due to differences in the input model GPDs so we neglect 
them). Henceforth, we suppress the factorized t-dependence and give predictions for t = 0. We will specify it for each 
case later. The LO and NLO coefficient functions, T^^^/^^ are taken from eqs. (14-17) of [Q and are summarized in 
appendix A of p^. They contain logarithms of the type log(l — X/C)"/(l — A^/C)"S with n, ni =0,1,2,3. Hence, 
depending on the region of integration, they can have both real and imaginary parts (i.e. if the argument of the log 
is positive or negative), which in turn generate real and imaginary parts of the DVCS amplitudes. 

In eq. (p3|) , we employ the +ie prescription through the Cauchy principal value prescription ("P.V.") which we 
implemented through the following algorithm: 

p.v. ^ dXT(^^-?j T{x, c, g2) = j^xrl^- 1^ {t{x, c, q^) - HQ, C, Q')) + 

dXTl^-^- 1^ {HX, C, 0') - ^(C, C, Q')) + ^(C, C, Q^) dXT 1^-^-1^. (44) 

Each term in eq. (|4j) is now either separately finite or only contains an integrable logarithmic singularity. This 
algorithm closely resembles the implementation of the -I- regularization in the evolution of PDFs and GPDs. We 
note that the first integral (in the ERBL region) is strictly real. The second and third terms contain both real and 
imaginary parts (which are generated in the DGLAP region). Explicit expressions for the real and imaginary parts 
of the DVCS amplitudes are given in Sec. (II) of [Q. They involve integrals over the coefficient functions which may 
be calculated explicitly and are given in Appendix B of [ p4[ . The real and imaginary parts of the unpolarized and 
polarized DVCS amplitudes were computed using a FORTRAN code based on numerical integration routines. We 
implemented the exact solution to the RNG equation for in LO or NLO in our calculation, as appropriate, to be 
consistent throughout our analysis. 



^^refering to vector and axial- vector currents for the quarks 
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B. Numerical results: unpolarized case 



Unpol. Re and Im singlet DVCS Amplitude in LO/NLO 
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FIG. 3. The Q^-dependence of the real and imaginary parts of the quark singlet DVCS amplitude. The solid (dashed) curve 
is the real part in LO (NLO) and the dotted (dashed-dotted) curve is the imaginary part in LO (NLO). 

In this subsection we illustrate the Q^- and ^-dependence, as well as the size, of the real and imaginary parts 
of the unpolarized DVCS amplitudes for the new input distributions (MRST99, CTEQ5M), calculated at LO and 
NLO accuracy. Fig. ^ shows the Q^-dependence of the real and imaginary parts of the quark singlet contribution, 
at LO and NLO, for two values oi — 0.1,0.0001, representative of HERMES and HERA kinematics, respectively. 
Correspondingly, Fig. ^ shows the gluon contributions, which start at NLO. Note the strong Q^-dependence in NLO 
of the imaginary part of the quark singlet and gluon amplitude at small C. This might raise concerns about the 
convergence of the perturbative expansion, especially when comparing NLO with LO in the quark singlet, where NLO 
grows much stronger with than LO. This is due to the same type of logarithmic divergences as X ^ ^ in both the 
evolution kernels and NLO coefficient functions. These divergences enhance the region around C, which is important 
for the value of the imaginary part psf , more quickly in than drops as increases. However, the quark 
singlet amplitude itself is not an observable quantity at NLO but rather the sum of quark singlet and gluon. When 
comparing the physical amplitudes at LO and NLO, the relative NLO corrections decrease as increases, as they 
indeed should psf . 

The above figures are complemented by Figs. H and Figs. ^, which show the C-dependence at fixed for the quark 
singlet and gluon contributions, respectively. Again we would like to point out the remarkable power-like behavior 
of the unpolarized amplitudes in C for fixed already remarked upon and explained in (with the exception of 
the NLO real parts of the quark singlet amplitude for CTEQ5M for the values plotted due to a somewhat strange 
combination of small quark input and comparatively large gluon input. At higher Q^, the CTEQ5M distribution also 
displays the characteristic power-like behavior in C^. 
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Unpol. Re and Im gluon DVCS Amplitude in NLO 
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FIG. 4. The Q'^-dependence of the real (solid line) and imaginary (dotted line) parts of the unpolarized gluon DVCS 
amplitude, which starts at NLO, for two representative values of (. 



Unpol. Re and Im singlet DVCS Amplitude in LO/NLO 
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FIG. 5. The real and imaginary parts of the unpolarized quark singlet DVCS amplitude, as a function of The solid 
(dashed) curve is the real part in LO (NLO) and the dotted (dashed-dotted) curve is the imaginary part in LO (NLO). To be 
able to plot the NLO real part of the quark singlet amplitude in a viewable manner, we removed the first point for CTEQ5M 
in the upper left plot (since it was negative and thus not easily handled in a log-log plot) and shifted the NLO curve of the real 
part in the upper right plot upward by an amount of C(l), again to avoid negative numbers. 
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Unpol. Re and Im gluon DVCS Amplitude in NLO 




FIG. 6. The real and imaginary parts of the unpolarized gluon DVCS amplitude as functions of (, for fixed Q . The solid 
curve is the modulus of the real part (which is actually negative) and the dotted curve is the modulus of the imaginary part of 
the gluon amplitude, which is also negative! 



C. Specification of t-dependence and lielicity flip amplitudes 

In this subsection we specify the ^-dependence of the various unpolarized/polarized helicity non- flip/flip DVCS 
amplitudes, Tii, Hi, £i, fi, for each parton species a — S{a = u,d,s),g. These are required in eqs.(H, ^, |^) to fully 
specify Tdvcs (our choices follow those of closely) 0. The various form factors are specified as follows: 

Q\t) = F^{t)T^'''{C Q^) , H?(C, Q^ t) = Gl{t)T^'\C Q2) , 
£:?(C, Q\t) = F^{t)T^;''{C, Q') , £^{C Q^ t) = G^r^^^iC, Q2) , (45) 

where, in the helicity flip case on the second line, we have introduced additional subscripts, £i,Ei, on the right hand 
side to distinguish this case from the helicity non-flip one considered explicitly above. For the up and down quark 
flavors we exploit the fact that proton and neutron form an iso-spin doublet to arrive at: 

2Fnt) = 2Ff{t) + Fr{t), 2Ffit) = F[it) + 2F'^{t) , for z = 1, 2 , (46) 

corresponding to the Dirac and Pauli form factors, respectively (see eqs.(p^, p^)). For the helicity non- flip polarized 
case we choose |^ : 

G^it) =(^-^) G-(t) = (1 - -L)-3 , (47) 

where niA = 0.9 GeV |^^. The strange quark and the gluon sea-like form factors were chosen in using the 
counting rules for clastic form factors which give for large t. Hence, for the unpolarized case the electric and 



In fact the GPDs themselves may already be considered to have a given helicity specification. The assumption of a factorized 
t-dependence for the GPDs allows the specification of the i-dependence to be moved to the amplitude level. 
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magnetic form factors were chosen to be 



1 + K'^'f 



-3 



(48) 



where we further assume k^^^ = 0. This gives (cf. eq.(|l2[)) 



-3 



0, 



(49) 



with my = 0.84 GeV. 

For the hehcity-flip case, we still have to specify the polarized and unpolarized GPDs. Unfortunately there are no 
inclusive analogs to guide us. However, we do know that the unpolarized helicity flip GPD contains a D-term with a 
relative minus sign as compared to the helicity non-flip GPD, such that when added they cancel |^ In contrast to |3|] 
we set the GPD equal to this D-term, which is non-zero only in the ERBL region^. We then evolve this distribution 
in and use it as an input to calculate T^^^ for use in eq.(|45|), with the for various a specified in eqs.(^, |9|). 

For the polarized helicity flip amplitude, 7^, it was observed that, in a similar fashion to the effective pseudo scalar 
form factors in /3-decay |33) , one can approximate it at small t by the pion pole (see for example |^ ) . Thus we chose 
the corresponding GPD to only contain the asymptotic pion distribution amplitude given by 



82~(X 



1 - 



X 



(50) 



in the ERBL region (and zero for X > Q. Since we use the asymptotic form we do not evolve the GPD and use it 



directly in the computation of the DVCS amplitude, T?''^[C,,Q'^). The i-dependence is given by the pion pole and 



thus we find 



^i(C,Q',i) 



^"^(C,Q'), i0Ya = u,d, 



(51) 



where g^'' = 1.267, M is the nucleon mass and T~'^ — 



£i 



for all Q^. For the s-quark and the gluon we set £ to 
zero. This completes the specification of the i-dependence of DVCS amplitudes, which arc then used to compute the 
various DVCS observables defined in Sec. ||. 

We close this section with a few comments. We note that £i,£i only have real parts since the GPD is zero for 
X > C,. Furthermore, in the asymptotic limit ^ oo for finite C, when all the partons have accumulated in the 
ERBL region, all amplitudes have zero imaginary part and a real part which is given by an asymptotic distribution 
amplitude in the ERBL region, thus the overall amplitude, and hence the appropriately scaled triple differential cross 
section, remain non-zero, in contrast to inclusive DIS. 



V. LO AND NLO RESULTS FOR DVCS OBSERVABLES 



In this section we present results for t he tr iple differential cross section and for various asymmetries (AAA, SSA, CA, 
UPLT, CADSFL) defined in subsection |II B| , in kinematics appropriate for the HI, ZEUS and HERMES experiments. 
We show the results as functions of t (for fixed C, Q^), of ( (for fixed t, Q^) and of (for fixed t, C) using the various 
input distributions defined in subsection p^II B . For our predictions in HERA kinematics we assume an e^P scattering 
with a proton energy of 920 GeV and a positron energy of 27.5 GeV (except of course for the charge asymmetries, 
which also use an electron probe of 27.5 GeV). 



■^'^This has to be the case since the n-th. moment in X of the sum is a polynomial of degree n-1 in ^, whereas the n-th. moment 
of the non-flip and flip GPDs separately are polynomials of degree n. The highest power of the polynomial in each case is 
generated by the D-term and thus they must cancel. 

Strictly speaking this violates the polynomiality condition since the helicity flip GPD is multiplied by the Pauli form factor 
rather than the Dirac form factor. However since the helicity flip term is numerically insignificant for DVCS observables this 
model is suflicient for phenomenological purposes. 
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A. The triple differential cross section 



In Fig. and Fig. [s] we show the triple differential cross section of eq.(^^, as a function of t at fixed C and Q^, for 
our four input distributions. We note that at the common point — A GeV^, GRV98, MRS A' and MRST99 are in 
close agreement, for small ( = x. In general CTEQ5M and MRST99 experience only moderate changes in going from 
LO to NLO. 

At small C = 2;, for GRV98 and MRSA' input distributions, the NLO corrections are much larger. As increases, 
at fixed Ci the spread of the predictions is seen to decrease. This is partly due to the evolution washing out the 
differences between the various input distributions, and partly due to an increased significance of the BH process. 
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FIG. 7. The triple differential cross section in t for fixed x — and . The solid (dotted) curve is the MRSA set in LO 
(NLO) and the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 8. The triple difFerential cross section in t for fixed x — and Q^. The solid (dotted) curve is the CTEQ5M set in LO 
(NLO) and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 




FIG. 9. The triple differential cross section in a; = for fixed t and . The solid (dotted) curve is the MRSA set in LO 
(NLO) and the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 



do/dxdQ'd|t| in nb/GeV' 
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FIG. 10. The triple differential cross section in a: = for fixed t and Q^. The solid (dotted) curve is the CTEQ5M set in LO 
(NLO) and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 

In Figs. |[ |l^, we employ logarithmic scales for HERA kinematics {x < 10~^) to illustrate the ^-behavior for fixed 
and t. It is interesting to note that we find the same type of power law behavior, C", for the triple differential 
cross section, which also includes the BH process, as we found previously for the unpolarized DVCS amplitudes (see 
Figs. I -|). Note that the kink for Q'^ = 9 GeV^ is due to going from a BH dominated region at x = 0.0001 (y > 0.8) 
to a region where DVCS dominates {x > 0.0005 and y < 0.2). 
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The fact that the GRV98 and MRS A' input scale is Qq_= 4 GeV^, imphes that the available Q^-range at small 
X is rather limited for these distributions. Hence, in Fig. |ll| we show only the Q^-dependence of the MRST99 and 
CTEQ5M sets, which start at Qq = 1 GeV^. For smaller Q^, where DVCS dominates BH, the cross section falls 
quickly with and one is sensitive to the details of the choice of input distribution. As increases and BH starts 
to dominate over DVCS, one observes that the curves begin to converge and one loses sensitivity to the details of the 
input GPD. 

We would like to point out that at large x all of the distributions produce fairly similar results and that the NLO 
corrections are tame. This is expected because each of the distributions have been strongly constrained by global fits 
to high statistics data in this region and hence behave similarly. Any observed differences may result in part from 
the holistic nature of the GPDs, which requires a continuous function for all X, both at the input scale and upon 
evolution. The real part of the amplitude is sensitive to an integral over the ERBL region X < C, so one may have 
some residual sensitivity to the behavior at very small X <C C: particularly if this behavior is extreme. The imaginary 
part is strongly influenced by the behavior a,i X — C, which in turn is constrained by a symmetry in the ERBL region. 
This sensitivity is further enhanced at smaller where the input PDFs are not yet so well constrained by direct 
(mainly inclusive) measurements. Hence, we optimistically speculate that a detailed measurement of the differential 
cross section, at HERA, for relatively low k, 1 — A GeV^, could help to discriminate between different standard 
input PDFs. 
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FIG. 11. The triple differential cross section in for fixed t and x — C^. The solid (dotted) curve is the CTEQ5M set in LO 
(NLO) and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 



B. The single spin asymmetries (SSA,UPLT) 

Having found a large spread of predictions for the triple differential cross section for different input GPDs, we now 
turn to two single spin asymmetries which are directly sensitive to the imaginary part of DVCS amplitudes. Let us 
first discuss the SSA, defined in cq. (pO[) , which can be measured using a polarized lepton probe on an unpolarized 
target. At small x the numerator in eq.([20|) is directly proportional to Im Hi (see the second term of eq.(30) of pl[). 

In Figs. |l^ we illustrate the t-dependence at fixed ^ and Q^, and note that the predictions range from as little 
as 5% to as much as 30%. We note that in general the t-dependence is rather flat for t > —0.1 GeV^, indicating that, 
to a certain extent, the i-dependence cancels in the ratio of eq.(^). This means that an experimental measurement of 
this asymmetry, even with a rather coarse binning in t, would be able to distinguish between different input scenarios, 
especially at larger values. At the common scale of = 4 GeV^, the GRV98 and MRSA' sets produce very 
similar numbers within LO and NLO, whereas CTEQ5M and MRST99 are more different within LO and NLO, at 
least at small x. The NLO to LO corrections are generally speaking small to moderate (5 — 30%). 



Figs. 14, ^ show that the SSA drops steeply in ^ = a: for flxed and t, suggesting that the HERA experiments 
will only be able to measure the SSA in the small ( regime £ [10^^, 10^^]). We would also like to point out that the 
differences for the different input sets becomes so small on the scale of the actual value of the SSA at larger ( in HERA 
kinematics, that only very high statistics would be able to discriminate between them. Thus, for HERA, only a small 
C measurement would gives a reasonable discrimination between the various inputs. For HERMES kinematics, where 
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for fixed ( = x,Q'^ one of course probes a different y = Q^/xS value, the SSA again becomes sizeable. Unfortunately, 
in this high ( = x regime, the SSA is proportional to a linear combination of imaginary parts of 7ii,?ii,fi rather 
than just Im Hi as at small ( (see second term of eq.(30) of pi|). 
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FIG. 12. The single spin asymmetry as a function of t, for fixed x — and Q^. The solid (dotted) curve is the MRSA set 
in LO (NLO) and the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 13. The SSA in t for fixed x and The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 
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Single Spin Asymmetry 
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FIG. 14. The SSA in x ^ ( for fixed t and Q^. The sohd (dotted) curve is the MRSA' set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 15. The SSA as a function of x = C for fixed t and Q^. The solid (dotted) curve is the CTEQ5M set in LO (NLO) and 
the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 

Finally, in Fig. |l^ we plot the Q^-dependence for fixed C — ^ For small ^ we observe that the magnitude 

of the SSA increases with Q'^, for both distributions. At large C, both MRST99 and CTEQ5M have very similar Q 
behavior and the NLO corrections are small. Note that at small ( the NLO corrections appear to grow in which at 
first sight looks strange. This is due to the fact that the BH process gains prominence relative to the DVCS process, so 
any differences between calculations of ImTLi for the interference term in the numerator become enhanced. However, 
the NLO corrections remain moderate between 15 — 30%. 
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Single Spin Asymmetry 
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FIG. 16. The SSA as a function of Q^ for fixed t und x = (. The solid (dotted) curve is the CTEQ5M set in LO (NLO) and 
the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 

If one switches now to a longitudinally polarized proton target as available at HERMES, but not currently planned 
for HERA, and an unpolarized lepton probe, one can form the unpolarized single spin asymmetry UPLT which is 
directly sensitive to the imaginary part of a combination of DVCS amplitudes. Furthermore, for the small regime 
within HERA kinematics, we find on inspection of the first term (proportional to A sine/)) of eq.(31) of |^ that the 
UPLT at small ^ is directly proportional to the imaginary part of the polarized amplitude Tii and that the other 
amplitudes like the numerically large Hi are suppressed by x. Hence, even though Im Hi is about a factor of one 
thousand smaller than Im H at x = 10^^, the suppression factor of x means that Im Hi constitutes only a 10% 
correction. Hence the UPLT is mainly sensitive to Im Hi at small x. Since this amplitude is numerically small, 
we found the UPLT asymmetry itself to be very small, at small x — f^nd thus virtually impossible to measure. 
Therefore, we do not show plots for HERA kinematics but rather only for HERMES kinematics. 
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FIG. 17. The UPLT as a function of t, for fixed x and Q^. The solid (dotted) curve is the MRSA' set in LO (NLO) and the 
dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 
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The UPLT Asymmetry 
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FIG. 18. The UPLT as a function of t, for fixed x and The sohd (dotted) curve is the CTEQ5M set in LO (NLO) and 
the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 

In Figs. Hand we show the UPLT as a function of t for fixed and ^ x = 0.15, i.e. in HERMES kinematics, 
where Im H and other ampHtudes contribute significantly. Hence, despite the fact that this asymmetry should be 
measurable at HERMES, the results would only be useful in the context of a program of asymmetry measurements 
that would enable the individual DVCS amplitudes to be isolated. Again a rather flat behavior is observed in t. We 
observe that all sets agree well with one another and their NLO corrections are small, in terms of percentages. In Figs. 
p^ , ^ we plot the UPLT as a function of C = a; at fixed t, Q^. Note the good agreement of all sets in the — A GeV'^ 
figures. This agreement does not bode well for the usefulness of the UPLT to discriminate between different input 
models. The NLO effects are found to be generally small. Finally, in Fig. |l] we show that the Q^-behavior of the 
UPLT asymmetry, is rather complicated (for the same reasons as in the case of the SSA in HERMES kinematics). 
Note the close agreement of all sets for the > 4 GeV^ behavior and the smallness of the NLO corrections. 
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FIG. 19. The UPLT as a function of a; = C, for fixed t and Q^. The solid (dotted) curve is the MRSA' set in LO (NLO) and 
the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 
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The UPLT Asymmetry 
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FIG. 20. The UPLT as a function oi x ^ (, for fixed t and Q^. The sohd (dotted) curve is the CTEQ5M set in LO (NLO) 
and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 
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FIG. 21. The UPLT as a function of for fixed t and x = (. In the left figure, the solid (dotted) curve is the CTEQ5M 
set in LO (NLO) and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO), and in the right figure, the solid 
(dotted) curve is the MRSA' set in LO (NLO) and the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 



C. The azimuthal angle asymmetry (AAA) 

We now discuss the (unpolarized) azimuthal angle asymmetry, defined in eq. (p9[), which directly probes the real 
part of DVCS amplitudes and only Re Tii at small C, — x where it dominates the other amplitudes (cf. the term 
proportional to cos0 in eq.(30) of |3^). Again we would like to point out that subtracting the BH contribution is 
necessary since it does not cancel in the asymmetry due to the 0-dependcncc of the propagators Vi and 1^2 ■ 

In general, Figs. |2^, reveal a rather flat behavior of the AAA in t, for fixed C,,Q^. We note that the spread 
of predictions coming from different inputs is rather large, indicating a strong sensitivity of the AAA to the input 
distributions and to the order in perturbation theory. The NLO corrections are very large, simply because the gluon 
enters at NLO for the first time, with a relative minus sign compared to the quarks in the real parts of the amplitudes. 
The wide spread in results indicates that the AAA is a highly sensitive discriminator between different input models, 
even those which have, up until now, agreed very well with one another. Hence, measuring the AAA both at HERA 
and HERMES with high precision is imperative for constraining the GPDs via a global fit. Note also that this strong 
sensitivity to the details of the shape and size of the GPD both in the ERBL and DGLAP region, as anticipated by 
the results in |2^,p3], indicates that an extraction of the GPDs with reasonable precision may be possible. 
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Azimuthal Angle Asymmetry 
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FIG. 22. The AAA in t for fixed x = C, and Q^. The sohd (dotted) curve is the MRSA set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 23. The AAA in t for fixed a; = C and The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 
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Azimuthal Angle Asymmetry 




-0.05 



-0.15 L 
10 
0.1 



0.05 


-0.05 



-0.1 



t = -0.25 GeV 



Q' = 4GeV' 



10"' 



t = -0.25 GeV 
Q' = 4GeV' 



-0.05 



-0.1 



-0.15 
10 10' 
0.1 



0.05 


-0.05 
-0.1 









/' 


: 


2 

t = -0.25 GeV 




Q' = 9GeV' 

1 


-4 


10"' 10 




1 ' 1 ' 

t = -0.25 GeV' 




- 6 25 GeV^ 











0.15 0.2 0.25 0.3 0.15 0.2 0.25 
X X 



0.3 



FIG. 24. The AAA in a; = C for fixed t and . The sohd (dotted) curve is the MRSA set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 25. The AAA in x = C for fixed t and . The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 

Figs. ^ and illustrate the strong decrease of the AAA as C, decreases within HERA kinematics (i.e. small 
for fixed ,t. However, as expected from dispersion relations, this behavior at small C of AAA, which is sensitive 
to ReJ-Li, is naturally not as steep as that of SSA, which is sensitive to ImTii- Note again the kink at = 9 GeV^ 
when going from a BH dominated region to a DVCS dominated one, i.e. from high to low y. 
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The ^-dependence of the AAA in the valence region probed at HERMES (where it has aheady measured the SSA) , 
shows very large NLO effects, consistent within all sets, which make the AAA large and thus measurable at HERMES. 
The large overall size and spread of predictions is encouraging for measurements of the AAA at large x since even 
here the discriminating power between GPD models is very good. 
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FIG. 26. The AAA in for fixed t and x = (. The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 

Finally we show the Q^-dependence of the AAA in Fig. (again only for MRST99 and CTEQ5M, due to their 
lower input scale). For small ( ~ x we see a fairly flat behavior of the AAA with with the NLO to LO changes 
staying fairly constant in Q^. At large C = x, the predictions for the two sets quickly approach one another in both 
LO and NLO as increases. Also the NLO to LO change decreases as increases in line with the argument that 
BH starts dominating at large due to the associated increase in y. 

D. The charge asymmetries (CA,CADSFL) 

Other asymmetries which measure the real part of DVCS amplitudes are the charge and charge double spin flip 
asymmetries (see eqs.(p3|, p^). These asymmetries require the measurement of DVCS with both the positron and the 
electron and, in the case of the CADSFL, a longitudinally polarized probe and target at the same time. Thus they 
are more difficult to measure than the AAA or other observables. The CA isolates the same combination of real parts 
of DVCS amplitudes as the AAA (i.e. the first term in eq.(30) of which is proportional to cos0), and would 
therefore serve as a useful complementary measurement. 

The CADSFL isolates the combination of real parts of DVCS amplitudes given in the second term of eq.(31) of 
The UPLT asymmetry isolates the imaginary part of the same linear combination. Thus a combined measurement of 
both CADSFL and UPLT at small x reveals information about the real and imaginary parts of the polarized helicity 
non-flip DVCS amplitude Hi- 

We discuss the CA first. Figs. p7[ - |3l| reveal that the CA is very similar to the AAA in the BH dominated region, 
and is within 20% of it in the DVCS dominated region for both LO and NLO. It can be easily seen that the t, x and 
behavior is very similar for both asymmetries. This behavior can be understood quite easily: in contrast to the 
AAA case, for CA the interference term drops out in the normalization in eq. ( p3| ) due to the sign change of the 
interference term in going from a positron to an electron probe. The numerator is the same for both asymmetries. 
Hence, deviations between the AAA and the CA results show directly the influence of the interference term relative 
to the BH and DVCS contributions. The larger the deviation the more important the interference term. Hence, a 
precise measurement of the CA and the AAA will serve as a very good consistency check of the experimental analysis, 
since strong deviations between the CA and AAA are not expected. One caveat here is the possible importance of 
higher twist contributions in the interference term. 
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FIG. 27. The CA in t for fixed x = C and Q^. Tlie solid (dotted) curve is the MRSA set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 28. The CA in t for fixed x = ( and The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 
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FIG. 29. The CA in x = C for fixed t and Q^. The sohd (dotted) curve is the MRSA' set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 
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FIG. 30. The CA in a; = C for fixed t and Q^. The sohd (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 
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FIG. 31. The CA in for fixed t and 2: = The soUd (dotted) curve is the CTEQ5M set in LO (NLO) and the dashed 
(dashed-dotted) curve is the MRST99 set in LO (NLO). 



Moving now to a longitudinally polarized probe and target, we will study the charge asymmetry in a double spin 
flip experiment, the CADSFL. 




FIG. 32. The CADSFL in t for fixed 2: = C and Q^. The sohd (dotted) curve is the MRSA' set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 

For HERA kinematics the CADSFL turns out to be too small to be measured, i.e. consistent with zero for any 
practical purposes, and thus we will not show plots for HERA kinematics. For HERMES kinematics, the asymmetry 
is of the order of a few percent, so one might hope to be able to measure it. In Figs. ^ and ^ we show the t- 
dependence for fixed ^ = x and Q^, which turns out to be rather steep. In Figs. ^ and ^ we show the ^-dependence 
for fixed i, . They indicate that although in percentage terms the spread of predictions is large the overall size of 
the CADSFL (only a few percent) would seem to indicate that exploiting this spread of predictions is impractical. 

Finally we show the Q^-behavior of CADSFL in Fig. for fixed t, C. We see quite a moderate to large difference 
in going from LO to NLO. Note the sign difference between the GRV98 and MRSA' sets on the one hand and the 
CTEQ5M and MRST99 sets on the other hand. This gives some hope to use this asymmetry as a discriminator 
between different model inputs. Generally speaking, due to its overall smallness, it is not as promising a candidate 
as the CA asymmetry as a good DVCS observable, even though it is directly sensitive to the real part of a polarized 
amplitude at small C. This concludes our presentation of DVCS observables in LO and NLO. 
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The CADSFL Asymmetry 
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FIG. 33. The CADSFL in t for fixed x ^ ( and Q^. The sohd (dotted) curve is the CTEQ5M set in LO (NLO) and the 
dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 
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FIG. 34. The CADSFL in x = C for fixed t and Q^. The sohd (dotted) curve is the MRSA' set in LO (NLO) and the dashed 
(dashed-dotted) curve is the GRV98 set in LO (NLO). 




FIG. 35. The CADSFL in x = C for fixed t and The solid (dotted) curve is the CTEQ5M set in LO (NLO) and the 
dashed (dashed-dotted) curve is the MRST99 set in LO (NLO). 
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FIG. 36. The CADSFL in for fixed t and x = (. In tfie left figure, tlie solid (dotted) curve is tlie CTEQ5M set in LO 
(NLO) and the dashed (dashed-dotted) curve is the MRST99 set in LO (NLO), and in the right figure, the solid (dotted) curve 
is the MRS A' set in LO (NLO) and the dashed (dashed-dotted) curve is the GRV98 set in LO (NLO). 



E. Summary of optimal DVCS observables 



We summarise the results of the previous subsections as follows: 

• The triple differential cross section is quite a good discriminator between different GPD models, particularly 
at values close to the input scale (at higher the BH process dominates). We predict a steep rise with 
decreasing reflecting the underlying powerlike behavior in ^ of the DVCS amplitudes. 

• The SSA is quite sizable and well behaved in NLO for all the input GPDs concerned, which makes it a very 
good candidate to be measured both at HERA and HERMES. For small ( (HERA kinematics) we predict a 
strong increase in magnitude as C decreases (a feature shared by AAA and CA) and good discriminating power 
between input models, especially at larger Q^. 

• The AAA and CA seem to be very good candidates for discriminating between different input GPDs simply 
because they measure the real part of DVCS amplitudes which are very sensitive to the underlying details of 
the GPDs, in Radyushkin's ansatz, especially at NLO where the gluon enters with a relative minus sign for the 
real part of its amplitude. Practical measurements for both asymmetries appear feasible for both HERMES and 
HERA kinematics. 

• The measurement of UPLT (CADSFL) appears to be practical only at large C in HERMES kinematics, where it 
tests a linear combination of imaginary (real) parts of several DVCS amplitudes. Therefore these measurements 
have very limited discriminating power between model GPDs. 

In summary it seems that four out of six observables are large enough to make measurements at HERA and HERMES 
feasible, and which have good discriminating power between input scenarios. Thus the prospects of unraveling the 
details of the DVCS process experimentally are in principle very good. 



VI. COMPARISON WITH EXPERIMENT AND OTHER CALCULATIONS 



In this section we compare our LO and NLO results for each of our input sets with the available published experi- 
mental data from HI and HERMES. The ZEUS collaboration announced the first measurement of DVCS at HERA 
(see also [Q), but has yet to publish a cross section. 
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HI recently published their data on the measured DVCS cross section on both the lepton level and the photon 
level. They use the equivalent photon approach, which relates the pure DVCS cross section on the lepton level (with 
the interference neglectedm and the pure BH term subtracted) to the virtual-photon proton cross section: 



dydQ^ =^-^vcsilP^lp) where r = . (52) 



By integrating eq. (^_7|) over t, 0, and changing variables from Xbj to y, we can establish the formula for the photon-level 
cross section in terms of our amplitudes as 

(^DVCs{l*P -> IP) = \TDVCS?\t=G, (53) 

where B stems from the i-integration and, within our model for the t-dependence, B ~ 6.5 GeV~^ (with a surprisingly 
small spread of about 1 GeV~^ !). We observed that for the kinematical region of the HI data which is limited to 
small X, i.e. to small C,, we could write 

dt ITcvcsP « hTDvcsf\t=o (54) 



as if we had assumed a global exponential dependence, e^*, as was used by HI in their comparison to other 
calculations (dropping all amplitudes except Ti). On numerical inspection we found that indeed Ti. clearly dominates 
all the other amplitudes for HERA kinematics. Therefore, in comparing to the data, we can safely neglect all other 
amplitudes (i.e. polarised and helicity-flip) in the DVCS square term. 

The same assumption was made in previous LO QCD jsj and two-component dipolc model calculations, 
which both reproduce the HI data quite well. Since ours is a QCD calculation we point out the differences and 
similarities with The latter was based on a LO input (CTEQ3L) with the assumption that the GPDs are equal to 
the PDFs at the input scale. The imaginary part of the unpolarized helicity non-flip amplitude was then computed 
using the aligned jet model and the imaginary part, rather than the GPD, was then evolved to higher Q^. This is 
equivalent to what we did at LO^ since the imaginary part of the DVCS amplitude is simply tt times the GPD at 
C — Xbj, legitimising the approach in ||]. The DVCS triple differential cross section, at small x, was then computed 
by comparing the imaginary part of DVCS to DIS and thus introducing the structure function F2 through the optical 
theorem. In order to reconstruct the real part of the DVCS amplitude a dispersion relation approach was used 
which exploited the slope of F2 in ln(l/a;) at small x (which was extracted from data to give rj = RcAdis/ IttiAdis) 
together with the comparative factor R — /to^dvcs / -^"^^dis • The specification was completed with the additional 
assumption of a global i-dependence, e^*. In contrast to we computed both parts of the DVCS amplitudes directly 
and used them in the unapproximated expressions for the DVCS triple differential cross section, which also contains 
the polarized as well as helicity flip amplitudes (however these can be safely neglected in the cross section at small 
C). Furthermore, we assumed a dipole type t-dependence which is close to the exponential behavior at small C and t, 
as well as being correct at larger t (and small C)- 

In Fig. |7|we plot aovcsil*?) = 4.5 GeV^ in W (this is of course similar to a plot in 1/x = W^/Q'^) for 

GRV and MRST input models at LO and NLO. Also shown is the HI data, for which they quote >= 4.5 GeV^, 
with systematic and statistical error bars added in quadrature. In Fig. we plot (Jjjycsi'l*?) in for fixed 
W = 75 GeV. The most obvious observation to make is that all of the the curves lie well above the data. However 
they do have the correct shapes (rising with W and falling rapidly with Q^). 



■^^This may be ignored at small y, typical of HERA kinematics, and after integration over 0. 

^''Note however that our LO calculation is based on NLO input PDFs, and a different ansatz for the GPDs, so we don't 
necessarily expect the results to agree closely. 
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FIG. 37. The photon level cross section a{-y*P 
at LO and NLO. Also shown are the recent HI 
quadrature. 



20 40 60 80 100 120 140 160 180 200 
W (GeV) 

7P) as a function of W at fixed 



4.5 GeV^ for GRV98 and MRST99 



data, at Q = 4.5 GeV , with systematic and statistical errors added in 



100 



O 
in 



10 



CO 



CO 

E 



LO GRV 
NLO GRV 
LO MRST 
NLO MRST 
HI data 



0.1 ^ ' ' ' ' ' ' ' 1 

2 4 6 8 10 12 14 16 

Q'^2 (GeV'^2) 

FIG. 38. The photon level cross section P -* 7P) as a function of at fixed W = 75 GeV, for GRV98 and MRST99 at 
LO and NLO. Also shown are the recent HI data, at = 75 GeV, with systematic and statistical errors added in quadrature. 

In order to understand this discrepancy we re-examined Radyushkin's input model in detail. For the unpolarized 
gluon GPD in |^ we found a rather moderate enhancement of about 20% of the GPD relative to the forward case 
at the point X = However, for the quark singlet, in which q rather than xq is used in the double distributions, 
we found this enhancement to be as large as a factor four! Such a large enhancement stems from the fact that the 
quark singlet distributions from the chosen input GPDs are very singular in the small x region {q{x) oc a;~^~^ ), 
and therefore lead to an integrable singularity in the double distribution at x' = (the shape function for the quark 
singlet is oc X for a; — > and this reduces the degree of the singularity to an integrable one). This leads to a large 
enhancement close to X — C, and a strong and unsatisfactory sensitivity to the extrapolation of the input PDFs to 
very small x, where they have not yet been measured (cf. eqs.( |35| , ^6| )). This is clearly a very unsatisfactory behaviour 
from a physical point of view and thus we are led to conclude that Radyushkin's input model should be used only 
with non-singular inputs. 

To illustrate the strong sensitivity of the results to the particular choice of input GPD we make a very simple change 
to the input model, namely we shift the argument of the PDFs in eq.(14) of |^l|] from x to x + This (admittedly 
rather 'ad hoc') modification, which we make purely for the purpose of demonstration, preserves the symmetries in 
the ERBL region and the forward limit, but can be expected to spoil the polynomial! ty properties, which are then 
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restored^ under evolution ||45[] . In making this modification we ensure that the PDFs are never sampled below x = 
(which removes the sensitivity to the very small x region) and can be expected to reduce considerably the enhancement 
of the GPDs at X = Figs. |3^ and |o| show the very dramatic effect on the cross section. The theory curves for the 
modified ansatz now undershoot the experimental data by about a factor of two and appear to reproduce the shape 
in W and rather well. In summary we claim that this illustrates that the HI data are already able to begin to 
constrain the input GPDs. It seems clear that if one is to use Radyushkin's ansatz at very small C the input PDFs 
should be non-singular. Alternatively, one must invent a new parameterization of the GPDs at the input scale that 
retains the required features (correct polynomiality and symmetry properties and faithful reproduction of the forward 
limit) without the problem of a strong sensitivity to the very small x region and associated large enhancement of the 
GPDs relative to the PDFs at the point X = (. This might be achieved by choosing a non-singular double distribution 
at a very low input scale, with the small X behavior generated by perturbative evolution. 
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FIG. 39. The photon level cross section aij'P ~* -yP) as a function of W at ^ 4.5 GeV^ for GRV98 and MRST99 at 
LO and NLO using the shifted ansatz. 
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FIG. 40. The photon level cross section a(^*P 
at LO and NLO, using the shifted ansatz. 
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in which the evolution is performed using moments, it is shown that expansion coefficients in an orthogonal basis 
allow only even powers of C, in the expansion, in line with the symmetry properties of GPDs. Similarly, in our numerical 
solution, which is carried out in x space rather than moment space, the non-polynomial pieces will be killed by the kernels 
under convolution (the kernels encode the GPD symmetries in their structure). 
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We now turn to the SSA as measured by the HERMES collaboration They define the SSA, weighted with 

sin (j>, as follows: 

SSA = 2r#sin0(d.T_rf,i) 

They quote two values: SSA = -0.18 ± 0.05 ± 0.05 and (SSA) = -0.23 ± 0.03 ± 0.04. The first, which assumes that 
their missing mass equals the proton mass, is quoted at the following average values: <x>= 0.11, <Q^>= 2.5 GeV^ 
and <t>= —0.27 GeV^. The second (average) value is the SSA integrated over the missing mass at the same average 
values of x, and t. At this point in x,Q'^,t for the HERMES definition we find SSA = -0.365 for CTEQ5M and 
SSA = -0.35 for MRST99. 

Given the fact that the input models used do not describe the data at small C, and that we find the same type 
of enhancement effect when comparing GPDs with forward PDFs at both small and large the failure to describe 
the HERMES data is not surprising. In fact one needs higher statistics over a wide kinematic range in several 
DVCS observables to be able to begin to tune the input GPDs using a fitting method. Only then can one start 
discriminating between different choices for the PDFs in input models. Furthermore, for HERMES data one needs to 
know the normalizations of numerator and the denominator producing the measured asymmetry to help to constrain 
the normalization of the individual amplitudes (in order to be able to make any sensible statements about a comparison 
between theory and experiment). 

In comparing with other calculations we show qualitative agreement with [^,^,^ and quantitative agreement with 
pT|,p[ wherever we used the same input distributions (MRS A' and GS(A) only). 



VII. CONCLUSIONS 



We have presented a detailed next-to-leading (NLO) QCD analysis of deeply virtual Compton scattering (DVCS) 
observables. We quantified the NLO corrections and established which observables have the best prospects to be 
measured accurately at HERA and HERMES (the triple differential cross section and the azimuthal angle (AAA), 
single spin (SSA) and charge (CA) asymmetries). We have demonstrated that such measurements would have good 
discriminating power between different scenarios for the generalized parton distributions (GPDs), by examining four 
cases. It turns out that AAA and CA, which test the real part of unpolarized DVCS amplitude at small x are the 
most sensitive to the choice of GPD model. 

We performed a comparison with presently available DVCS data and showed that within Radyushkin's ansatz for 
the GPDs, all of the models we examined overshoot the HI data. This is due to their singular nature of the quark 
singlet at small x which leads to a strong sensitivity to the PDFs in the extremely small x region, where they have not 
yet been measured. We illustrated this strong sensitivity by artificially shifting the argument of the PDFs by C, which 
led to a large reduction in the theoretical predictions for the photon-level cross section (which then undershoots the 
data). This illustrates an urgent need for improved input models, which do not rely on a singular double distribution 
and also for more, high precision data to help constrain them. 
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